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Abstract. In this talk I review recent progresses in the detection of scalar 
. ^ I gravitational waves. Furthermore, in the framework of the Jordan-Brans- 

Dicke theory, I compute the signal to noise ratio for a resonant mass detector 
G^ ■ of spherical shape and for binary sources and collapsing stars. Finally I 

compare these results with those obtained from laser interferometers and 
from Einsteinian gravity. 
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The efforts aimed at the detection of gravitational waves (GW) started more than a 
quarter of century ago and have been, so far, unsuccessful [jl], Resonant bars have 
bJ). proved their reliability, being capable of continous data gathering for long periods of 

time 0, ^. Their energy sensitivity has improved of more than four orders of magnitude 
since Weber's pioneering experiment. But a further improvement is still necessary to 



X 



}-i ' achieve successful detection. While further developments of bar detectors are under way, 

^ two new generations of earth based experiments have been proposed: detectors based 

on large laser interferometers are already under construction 0, resonant detectors of 
spherical shape are under study [0. 

In this lecture I report on a series of papers in which the opportunity of intro- 
ducing resonant mass detectors of spherical shape was studied. As a general motivation 
for their study, spherical detectors have the advantage over bar-shaped detectors of a 
larger degree of symmetry. This translates into the possibility of building detectors of 
greater mass and consequently of higher cross section. 

Besides this obvious observation, the higher degree of symmetry enjoyed by the 
spherical shape puts such a detector in the unique position of being able of detecting 
GW's with a spin content different from 2. This means to test non-Einsteinian theories 
of gravity. 

I would in fact now like to remind the reader of the very special position of Ein- 
stein's general relativity (GR) among the possible gravitational theories. Theories of 
gravitation, in fact, can be divided into two families: metric and non-metric theories 0. 
The former can be defined to be all theories obeying the following three postulates: 

^ E-mail: Fucito@roma2.infn.it. 



1 



• spacetime is endowed with a metric; 

• the world hnes of test particles are geodesic of the above mentioned metric; 

• in local free-falling frames, the non-gravitational laws of physics are those of special 



It is an obvious consequence of these postulates that a metric theory obeys the principle 
of equivalence. More succintly a theory is said to be metric if the action of gravitation on 
the matter sector is due exclusively to the metric tensor. GR is the most famous example 
of a metric theory. Kaluza-Klein type theories, also belong to this class along with the 
Brans-Dicke theory. Different representatives of this class differ for their equations of 
motion which in turn can be deduced from a lagrangian principle. Since there seems 
to be no compelling experimental or theoretical reasons to introduce non-Einsteinian 
or non-metric theories, they are sometimes considered a curiosity. This point should 
perhaps be reconsidered. String theories are in fact the most serious candidate for a 
theory of quantum gravity, the standard cosmological model has been emended with 
the introduction of inflation and even the introduction of a cosmological constant (which 
seems to be needed to explain recent cosmological data) could imply the existence of other 
gravitationally coupled fields. In all of the above cited cases I am forced to introduce 
fields which are non-metrically coupled in the sense explained above. 

In the first section of this lecture I will explain that a spherical detector is able to 
detect any spin component of an impinging GW. Moreover its vibrational eigenvalues 
can be divided into two sets called spheroidal and toroidal. Only the first set couples 
to the metric. This leads to the opportunity of using such a detector as a veto for non- 
Einsteinian theories. In the second section I take as a model the Jordan-Brans-Dicke 
(JBD), in which along with the metric I also have a scalar field which is metrically 
coupled. I am then able to study the signal to noise ratio for sources such as binary 
systems and collapsing stars and compare the strenght of the scalar signal with respect 
to the tensor one. Finally in the third and last section I repeat this computation in the 
case of the hollow sphere which seems to be the detector which is most likely to be built. 

2. Testing Theories of Gravity 

2.1. Free Vibrations of an Elastic Sphere 

Before discussing the interaction with an external GW field, let us consider the basic 
equations governing the free vibrations of a perfectly homogeneous, isotropic sphere of 
radius R, made of a material having density p and Lame coefficients A and n [§]. 

Following the notation of 0, let Xi,i = 1,2,3 be the equilibrium position of the 
element of the elastic sphere and x- be the deformed position then the 
displacement vector. Such vector is assumed small, so that the linear theory of elasticity 
is applicable. The strain tensor is defined as Uij = (1/2) (wj ,,- + Uj^i) and is related to 
the stress tensor by (Jij = 6ijXuii + 2fiUij. The equations of motion of the free vibrating 
sphere are thus 



relativity. 




(2.1) 
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with the boundary condition: 

n^aij = (2.2) 

dX r = R where rii = Xi/r is the unit normal. These conditions simply state that the 
surface of the sphere is free to vibrate. The displacement Ui is a time-dependent vector, 
whose time dependence can be factorised as Ui{x,t) = Ui{x)exp{iLjt), where u is the 
frequency. The equations of motion then become: 

fiV'^Ui + (A + iiy jUj) = —co'^pui (2.3) 

Their solutions can be expressed as a sum of a longitudinal and two transverse vectors 



13: 

uix) = CoV0(f) + CiLxix) + C2V X Lx{x) (2.4) 

where Co,Ci,C2 are dimensioned constants and L = x x V is the angular momentum 
operator. Regularity at r = restricts the scalar functions and x to be expressed as 
(f){r,9,(p) = ji{qr)Yim{9,ip) and x{'r;0,(p) = ji{kr)Yim{9,ip). Yim{9,ip) are the spherical 
harmonics and ji the spherical Bessel functions ||1 1|| : 



/ 1 \ ' / sinx\ 

= pc<j^/(A + 2/i) and k"^ = puj'^/fi are the longitudinal and transverse wave vectors 
respectively. 

Imposing the boundary conditions (|2.2|) at r = i? yields two families of solutions: 



Toroidal modes: these are obtained by setting Cq = C2 = 0, and Ci 7^ 0. In this 
case the displacements in (|2.4|) can be written in terms of the basis: 



iZimir,e,^)=TUr)LYiUO,^) (2.6) 

with Tniir) proportional to ji{knir). The eigenfrequencies are determined by the 
boundary conditions (|2.2| ) which read 



h{kR) = Q (2.7) 

where 

(2.8) 



Spheroidal modes: these are obtained by setting Ci = 0, Co 7^ and C2 7^ 0. The 
displacements of (p.4|) can be expanded in the basis 



^^r.im{x) = Ani{r)Yi^{e, ip)n - Bniir)n x LYi„(^, cp) (2.9) 

where Ani{r) and Bni{r) are dimensionless radial eigenfunctions which can be 
expressed in terms of the spherical Bessel functions and their derivatives. The 
eigenfrequencies are determined by the boundary conditions ( |2.2| ) which read [^] 

where 2 



The eigenfrequencies can be determined numerically for both toroidal and spheroidal 
vibrations. Each mode of order / is {21 + l)-fold degenerate. The eigenfrequency values 
can be obtained from : 

[JI{kR)ni 



p R 

2.2. Interaction of a Metric GW with the Sphere Vibrational Modes 



(2.12) 



The detector is assumed to be non-relativistic (with sound velocity f ^ <^ c and radius 
R <^ \ the GW wavelength) and endowed with a high quality factor {Qni = oJniTni ^ 1, 
where Tni is the decay time of the mode nl). The displacement -u of a point in the detector 
can be decomposed in normal modes as: 



U{x,t) = ^AN{t)tljN{x) 

N 



(2.13) 



where collectively denotes the set of quantum numbers identifying the mode. The 
basic equation governing the response of the detector is 



^N{t) + r^'i7v(t) + ujj,AN{t) = fN{t) 



(2.14) 



I assume that the gravitational interaction obeys the principle of equivalence which 
has been experimentally supported to high accuracy. In terms of the so-called electric 
components of the Riemann tensor Eij = Roioj, the driving force fN{t) is then given by 



m 



f^it) = -M-^E,,{t) / i;'^{x)x^ pd'x 



(2.15) 



where M is the sphere mass and I consider the density p as a constant. In any metric 
theory of gravity Eij is a 3 x 3 symmetric tensor, which depends on time, but not on 
spatial coordinates. 

Let us now investigate which sphere eigenmodes can be excited by a metric GW, 
i.e. which sets of quantum numbers N give a non-zero driving force, 
a) Toroidal modes 

The eigenmode vector, ipnim can be expressed as in eq. (|2.6| ). Up to an adimensional 
normalisation constant C, the driving force is 
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Using the equations 
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(2.17) 



(2.18) 



A-n (/ + m) 

the integration over can be performed. Eq. ( p.l6|) then contains integrals over 9 of the 



form: 
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(2.20) 



After integration by parts, the derivative terms in eqs. (|2.19| ) and ( |2.20|) exactly cancel 
the non-derivative ones. The remaining boundary terms vanish too, thanks to the peri- 
odicity of the trigonometric functions and to the regularity of the associated Legendre 
polynomials. The vanishing of the above integrals has a profound physical consequence. 
It means that in any metric theory of gravity the toroidal modes of the sphere cannot 
be excited by GW and can thus be used as a veto in the detection. 

b) Spheroidal modes 
The forcing term is given by: 



fi\t) = -M-'E,,{t) / xH-AN{r)Yi^{9,v) - BN{r)e'^'—LuYi^{9,^) ]pd'x (2.21) 



X 



ink -^n 



One is thus lead to compute integrals of the following types 

j x^x'Yi„,{9,ip)d^x 

and 

x^x'LkYi^{9,ip)d^x 



(2.22) 



(2.23) 



Since the product x^x^ can be expressed in terms of the spherical harmonics with I = 0,2 
and the angular momentum operator does not change the value of /, one immediately 
concludes that in any metric theory of gravity only the / = 0, 2 spheroidal modes of the 
sphere can be excited. At the lowest level there are a total of five plus one independent 
spheroidal modes that can be used for GW detection and study. 

2.3. Measurements of the Sphere Vibrations and Wave Polarization States 

From the analysis of the spheroidal modes active for metric GW, I now want to infer 
the field content of the theory. For this purpose it is convenient to express the Riemann 
tensor in a null (Newman- Penrose) tetrad basis [0. 

To lowest non-trivial order in the perturbation the six indipendent "electric" com- 
ponents of the Riemann tensor may be expressed in terms of the Newmann- Penrose (NP) 
parameters as 

-i?e^4-$22 /m^4 -2^/2Re^3' 
Eii = I /m^4 Re^4 - $22 2V2/m^3 ] (2.24) 
-2V2Re^3 2y2/m*3 -6^2 
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The NP parameters allow the identification of the spin content of the metric theory 
responsible for the generation of the wave 0. The classification can be summarized in 
order of incresing complexity as follows: 

• General Relativity (spin 2): \E'4 7^ while \1>2 = vE's = $22 = 0. 

• Tensor-scalar theories (spin 2 and 0): \E'4 7^ 0, \E'3 = 0, \1'2 7^ and/or $22 7^ {e.g. 
Brans-Dicke theory with ^1^4 7^ 0, \I^2 = 0, \I^3 = and $22 7^ 0). 

• Tensor- vector theories (spin 2 and 1): \E'4 7^ 0, \E'3 7^ 0, $22 = ^^2 = 0. 

• Most General Metric Theory (spin 2, 1 and 0): ^47^ 0, ^2 7^ 0, ^3 7^ and 
$22 7^ 0, {e.g. Kaluza-Klein theories with \E'4 7^ 0, \I^3 7^ 0, $22 7^ while \I^2 = 0). 

In eq. ( |2.24| ), I have assumed that the wave comes from a localized source with 
wave vector k parallel to the z axis of the detector frame. In this case the NP parameters 
(and thus the wave polarisation states) can be uniquely determined by monitoring the 
six lowest spheroidal modes. If the direction of the incoming wave is not known two more 
unknowns appear in the problem, i.e. the two angles of rotation of the detector frame 
needed to align k along the z axis. In order to dispose of this problem one can envisage 
the possibility of combining the pieces of information from an array of detectors \T^\. I 
restrict my attention to the simplest case in which the source direction is known. 

In order to infer the value of the NP parameters from the measurements of the 
excited vibrational modes of the sphere, I decompose Eij in terms of spherical harmonics. 
In fact, the experimental measurements give the vibrational amplitudes of the sphere 
modes which are also naturally expanded in the above basis. The use of the same 
basis makes the connection between the NP parameters and the measured amplitudes 
straightforward. In formulae 

E^At) = j:^Ut)S^j'^^ (2.25) 

l,m 

where S^j'^^ = 6ij/y/4n (with 6ij the Kronecker symbol) and S^f"^^ {m = —2, ..2) are five 
linearly independent symmetric and traceless matrices such as 

5g'-)nV = y;^, 1 = 0,2 (2.26) 

The vector rii in eqs. ( |2.26|) has been defined after eq. ( |2.2|) . 
Taking the scalar product I find 



C0,0(t) = Y5°/Ey(t) 

^TT „2 
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and then for the NP parameters 



C2,Ut) = T^StpE,,(t) (2.27) 



$22 = \/l|^^2.o(t)-y^Co,o(t) *2 = -^^^C2,o(t)-^^ico,o(t) 
/ X5 / X5 
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1/15 i 15 

Re^3 = T^W — [C2.i -C2,-i] /m^3 = T^W — [C2,i + C2,-i] (2.28) 
lb V 71" ' lb V 

Eqs. ( p.28| ) relate the measurable quantities with the GW polarization states, de- 
scribed by the NP parameters. Eq. ( |2.28| ) can be put in correspondence with the output 
of experimental measurements if the are substituted with their Fourier components 
at the quadrupole and monopole resonant frequencies which, for the sake of simplicity, 
I collectively denote by uq. The Q,m(i^o) can be determined in the following way: once 
the Fourier amplitudes A]^{ujo) are measured, by Fourier transforming ( p. 141) and ( |2.15|) 



I get the Riemann amphtudes Eij{uJo) which, using (|2.27|) , yield the desired result. 

In order to determine the An{ujo) amplitudes from a given GW signal the following 
two conditions must be fulfilled: 

• the vibrational states of the five-fold degenerate quadrupole and monopole modes 
must be determined. The quadrupole modes can be studied by properly combining 
the outputs of a set of at least five motion sensors placed in independent positions 
on the sphere surface. Explicit formulas for practical and elegant configurations 



of the motion sensors have been reported by various authors [|I5| |Ib[. The vibra- 
tional state of the monopole mode is provided directly by the output of any of the 
above mentioned motion sensors. If resonant motion sensors are used, since the 
quadrupole and monopole states resonate at different frequencies, a sixth sensor is 
needed. 

The spectrum of the GW signal must be sufficiently broadband to overlap with the 
antenna quadrupole and monopole frequencies. 



3. Gravitational Wave Radiation in the Jordan-Brans-Dicke Theory 

In this section I analyze the signal emitted by a compact binary system in the Jordan- 
Brans-Dicke theory. I compute the scalar and tensor components of the power radiated 
by the source and study the scalar waveform. Eventually I consider the detectability of 
the scalar component of the radiation by interferometers and resonant-mass detectors. 



3.1. Scalar and Tensor GWs in the JBD Theory 



In the Jordan-Fierz frame, in which the scalar field mixes with the metric but decouples 
from matter, the action reads |17| 



S 



^grav 
IbTT 



BD_ ILV 



+ ^ J d'^xL^[lpra,gf,u] , (3.1) 



where cu^^ is a dimensionless constant, whose lower bound is fixed to be u^^ ~ bOO by 
experimental data [|1^], g^^i, is the metric tensor, is a scalar field, and ip^a collectively 
denotes the matter fields of the theory. 
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As a preliminary analysis, I perform a weak field approximation around the back- 
ground given by a Minkowskian metric and a constant expectation value for the scalar 
field 

ip = ipo + ^ . (3.2) 



The standard parametrization ipo = 2{ujg^ + 2)/G{2lj^^ + 3), with G the Newton con- 
stant, reproduces GR in the limit a;^^ oo, which implies (po ^ 1/G. Defining the new 
field J 

dfMu = h^^- -rj^^h-r]^^ — (3.3) 

where h is the trace of the fluctuation h^^^, and choosing the gauge 

d^e"" = (3.4) 

one can write the field equations in the following form 

= -—V. (3.5) 

= , , S (3.6) 



2u;,,+3 



where 



^-(T^. + V) (3.7) 



2(2^^^ + 3)1 2 ifo 16n 



2 ipo 



(3.8) 



In the equation (|377| ), T^,y is the matter stress-energy tensor and is the gravitational 
stress-energy pseudo-tensor, that is a function of quadratic order in the weak gravita- 
tional fields 6fj,iy and ^. The reason why I have written the field equations at the quadratic 
order in 6*^;^ and ^ is that in this way, as I will see later, the expressions for 9^i, and ^ 
include all the terms of order (f /c)^, where v is the typical velocity of the source (New- 
tonian approximation). 

Let us now compute and S at the order (f/c)^. Introducing the Newtonian 
potential U produced by the rest-mass density p 

U{x,t)= [ -P^23_£x' (3.9) 

J \ X — X \ 

the total pressure p and the specific energy density 11 (that is the ratio of energy density 
to rest-mass density) I get (for a more detailed derivation, see P): 

= -p , (3.10) 

T (I t\ 

2(2cu^^ + 3)V 2 ip,) 

P ^l + n-3g + '^- + V') (3.11) 



2(2^,, + 3 V P ^Bz, + 2 
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Far from the source, the equations ( p.5|) and ( |3.6| ) admit wave-hke solutions, which 
are superpositions of terms of the form 



= A^^{x, uj) exp{ik°'Xa) + c.c. 



(3.12) 
(3.13) 



Without affecting the gauge condition (^.41), one can impose h = —2^/ipo (so that 9^ 



h^jj). Gauging away the superflous components, one can write the amphtude A^j, in 
terms of the three degrees of freedom corrisponding to states with hehcities ±2 and 



19|. For a wave traveUing in the z-direction, one thus obtains 
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(3.14) 



where b = B/lpq. 

3.2. Power emitted in GWs 

The power emitted by a source in GWs depends on the stress-energy pseudo-tensor t^^ 
according to the following expression 



P^^ = r'^ j MVl = r^ J <t^^ > Xk dVt 



(3.15) 



where r is the radius of a sphere which contains the source, VL is the solid angle, $ is the 
energy flux and the symbol < ... > implies an average over a region of size much larger 
than the wavelength of the GW. At the quadratic order in the weak fields I find 



< to, > 



.foc_ 
' 32n 



^^""""^"^ < (doOidoO > + < {doKp){doh^P) > 



(3.16) 



Substituting (|1|), ( gH into (|3l6|) , one gets 

< to, >= 
and using ( |3.14|) 



4 2 

IQtt 



Vo 



(3.17) 
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^(poc UJ 
' Stt 



en P + I ei2 P +(2cj^^ + 3) | b 



(3.18) 



From ( p.l8| ) I see that the purely scalar contribution, associated to 6, and the trace- 
less tensorial contribution, associated to e^^, are completely decoupled and can thus be 
treated independently. 

3.3. Power emitted in scalar GWs 

I now rewrite the scalar wave solution (|3.13| ) in the following way 



^{x,t) = ^{x,u)e-''^' + C.C. 



(3.19) 
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In vacuo, the spatial part of the previous solution ( |3.19|) satisfies the Helmholtz equation 

(V2 + cj2)^(f,cj) = (3.20) 
The solution of ( |3.20| ) can be written as 

00) = Y.X,mhf\ujr)Y,UO,^) (3.21) 

jm 

where h^j^\x) are the spherical Hankel functions of the first kind, r is the distance of 
the source from the observer, Yjjn{0, (p) are the scalar spherical harmonics and the coeffi- 
cients Xjm give the amplitudes of the various multipoles which are present in the scalar 
radiation field. Solving the inhomogeneous wave equation (^.6|), I find 



Xjm = IQttiuj / jiiujr')Y;„Xd,v)Six,uj) dV (3.22) 



V 

where ji{x) are the spherical Bessel functions and r' is a radial coordinate which assumes 
its values in the volume V occupied by the source. 

Substituting ( p.l6|) in (p.l5|) , considering the expressions (|3.19|) and ( |3.21|) , and 
averaging over time, one finally obtains 

Pscal ~ 7^ ^ ^, I -^jm I (3.23) 

To compute the power radiated in scalar GWs, one has to determine the coefficients 
Xjm, defined in ( p.22|) . The detailed calculations can be found in the appendix A of the 
third reference in |p, while here I only give the final results. Introducing the reduced 
mass of the binary system /i = 17117712/171 and the gravitational self-energy for the body 
a (with a = 1,2) 

n. = --[ #^ d'xd'x' (3.24) 

2 JVa \ X — x' \ 

one can write the Fourier components with frequency nuo in the Newtonian approxima- 
tion 



16V27r tujQipo 

(Xoo)n = — nJn{ne) (3.25) 

3 uo +2 a 
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T 2(1 - e2)V2[(i _ e')Une)/e' - J'Sne) / {ne)]} 



Substituting these expressions in (|3.23|) , leads to the power radiated in scalar 
the n-th harmonic 

where the monopole, dipole and quadrupole terms are respectively 



(3.28) 
GWs in 
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15(cj^^ + 2) a^c^ 



g{n;e) 



terms 



(3.31) 
(3.32) 

The total power radiated in scalar GWs by a binary system is the sum of three 

(3.33) 



where 



pj=0 

pj=2 
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Pscal = P'=' + P'=' + P'=^ 

G4 2 2 2 
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Note that P^^^, P^^^, P-'=2 gQ ^ero in the limit cj^^ 



(3.34) 
(3.35) 
(3.36) 



00. 
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3.4- Scalar GWs 



I now give the exphcit form of the scalar GWs radiated by a binary system. To this end, 
note that the major semi-axis, a, is related to the total energy, E, of the system through 
the following equation 

a = (3.37) 

Let us consider the case of a circular orbit, remembering that in the last phase of evolution 
of a binary system this condition is usually satisfied. Furthermore I will also assume 
mi = 1712- With these positions only the quadrupole term, (|3.32|) , of the gravitational 
radiation is different from zero. The total power radiated in GWs, averaged over time, 
is then given by ( |3.34|) -( ^736|) 



P 



+ 3) + 1] 



(3.38) 



15 (cj^^ + 2) c5 rfs 

where d is the relative distance between the two stars. The time variation of d in one 
orbital period is 



d 



Gmim2 
2E^ ' 



(3.39) 



Finally, substituting ( |3.37| ), ( |3.38| ) in ( |3.39| ) and integrating over time, one obtains 

(1 12cu„„ + igC^mimzmX^^^ , 



d 



,15 ^sz, + 2 



(3.40) 



where I have defined t = tc — t, tc being the time of the collapse between the two bodies. 

From ( p.21|) , ( p.25|) -( pl28D and one can deduce the form of the scalar field (see 
appendix B of the third in for details) which, for equal masses, is 

V H — - — [h ■ v) 
d 



m 



(3.41) 



,2 I (c. 

where r is the distance of the source from the observer, and h is the versor of the 
line of sight from the observer to the binary system center of mass. Indicating with 
7 the inclination angle, that is the angle between the orbital plane and the reference 
plane (defined to be a plane perpendicular to the line of sight), and with the true 
anomaly, that is the angle between d and the x-axis in the orbital plane x-y, yields 
n ■ d = c? sin 7 sin 'i/'. Then from ( |3.41| ) one obtains 



which can also be written as 



20(1171 



{2uj^^ + 3)c^dr 



sin^ 7 cos(2'?/'(t)) 



e(r) = eo(r) sin(x(r) + x) 
where x is an arbitrary phase and the amplitude ^o(''") is given by 

2Gfim 



(3.42) 



(3.43) 



{2uj^^ + Sydr 



sin^ 7 



1/4 
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V2cii 



1/4 



sm 7 



2{2uj^^+3)r Vl2^,, 
In the last expression, I have introduced the definition of the chirp mass M, 



(3.44) 
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3.5. Detectability of the scalar GWs 



Let me now study the interaction of the scalar GWs a spherical GW detector. 

As usual, I characterize the sensitivity of the detector by the spectral density of 
strain Sh{f) [Hz]^^. The optimum performance of a detector is obtained by filtering the 
output with a filter matched to the signal. The energy signal-to-noise ratio SNR of the 
filter output is given by the well-known formula: 

r+oo \H(f)\'^ 

SNR = / df (3.45) 

where H{f) is the Fourier transform of the scalar gravitational waveform hs{t) = G^Q{t). 

I must now take into account the astrophysical restrictions on the validity of the 
waveform ( p.43|) which is obtained in the Newtonian approximation for point-like masses. 



In the following, I will take the point of view that this approximation breaks down when 



there are five cycles remaining to collapse |20, 21 



The five-cycles limit will be used to restrict the range of M^. over which my analysis 
will be performed. From ( |3.40| ), one can obtain 




Integrating (|3.46| ) yields the amount of phase until coalescence 

Xir) ^ f + ^ (3 47) 



Setting ( p.47| ) equal to the limit period, T5 cycles = 5(27r), solving for r and using ( p.46| ) 
leads to 

= 2.(6870 Hz) f^i^^)'^|^ (3.48) 



Taking = 600, the previous limit reads 



cycles = 27t{15A7Rz)^ (3.49) 

A GW excites those vibrational modes of a resonant body having the proper simmetry. 
In the framework of the JBD theory the spheroidal modes with / = 2 and / = are 
sensitive to the incoming GW. Thanks to its multimode nature, a single sphere is capable 
of detecting GW's from all directions and polarizations. I now evaluate the SNR of a 
resonant-mass detector of spherical shape for its quadrupole mode with m = and its 
monopole mode. In a resonant-mass detector, Sh{f) is a resonant curve and can be 
characterized by its value at resonance Sh{fn) and by its half height width Sh{fn) 
can thus be written as 

Skifn) = (3.50) 

Here (T„ is the cross-section associated with the n-th resonant mode, T is the ther- 
modynamic temperature of the detector and Qn is the quality factor of the mode. 
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The half height width of Sh{f) gives the bandwidth of the resonant mode 



A/„ = ^T-'/' (3.51) 

Here, r„ is the ratio of the wideband noise in the n-th resonance bandwidth to the 
narrowband noise. 

From the resonant-mass detector viewpoint, the chirp signal can be treated as 
a transient GW, depositing energy in a time-scale short with respect to the detector 
damping time. I can then consider constant the Fourier transform of the waveform 
within the band of the detector and write [^] 

^h[jn) 

The cross-sections associated to the vibrational modes with / = and / = 2, m = 
are respectively p 

a(„o) = ^n,3(^^^+2) ^^-^^^ 

Fn GMVg' 

<y{n2) = ——, 3.54 

6 c^{uJsD + 2) 

All parameters entering the previous equation refer to the detector M is its mass, Vs the 
sound velocity and the constants Hn and F„ are given in 0. The signal-to noise ratio 
can be calculated analytically by approximating the waveform with a truncated Taylor 
expansion around t = 0, where ojgit = 0) = Uni fl^, pO | 

1 ( duj 



hs{t) ^ Gio{t = 0) sin 



'^"''+2 lift I 



t=0 J 



(3.55) 



Using quantum limited readout systems, one finally obtains 

5 ■ 



{SNRn, 



1=0 



3 



32{u^^+2){12iu^^ + 19)hc 

sm 7 (3.5oj 



{SNR 



5 ■ 21/3F 



n 1=2 



— iT^sin 7 3.57 

which are respectively the signal-to-noise ratio for the modes with I = and / = 2, m = 
of a spherical detector. 

It has been proposed to realize spherical detectors with 3 meters diameter, made 
of copper alloys, with mass of the order of 100 tons This proposed detector has 

resonant frequencies of uji2 = 2tt ■ 807 rad/s and uio = 2tt ■ 1655 rad/s. In the case of 
optimally oriented orbits (inclination angle 7 = vr/2) and cj^^ = 600, the inspiralling of 
two compact objects of 1.4 solar masses each will then be detected with SNR = 1 up to 
a source distance r(u;io) — 30 kpc and r{ijji2) — 30 kpc. 
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4. The hollow sphere 



An appealing variant of the massive sphere is a hollow sphere [^. The latter has the 
remarkable property that it enables the detector to monitor GW signals in a significantly 
lower frequency range — down to about 200 Hz — than its massive counterpart for com- 
parable sphere masses. This can be considered a positive advantage for a future world 
wide network of GW detectors, as the sensitivity range of such antenna overlaps with 
that of the large scale interferometers, now in a rather advanced state of construction 
p5| , p6| . In this Section I study the response of such a detector to the GW energy emitted 
by a binary system constituted of stars of masses of the order of the solar mass. A hollow 
sphere obviously has the same symmetry of the massive one, so the general structure of 
its normal modes of vibration is very similar to that of the solid sphere. In particular, 
the hollow sphere is very well adapted to sense and monitor the presence of scalar modes 
in the incoming GW signal. The extension of the analysis of the previous Sections to 
a hollow sphere is quite straightforward and in the following I will only give the main 
results. Due to the different geometry, the vibrational modes of a hollow sphere differ 
from those studied in Section 2. In the case of a hollow sphere, I have two boundaries 
given by the outer and the inner surfaces of the solid itself. I use the notation a for the 
inner radius, and R for the outer radius. The boundary conditions are thus expressed by 



(JijUj 







at r = R and at 



{R>a> 0), 



(4.58) 



(|2.3| ) must now be solved subject to this boundary conditions. The solution that leads to 
spheroidal modes is still ( |2.9| ) where the radial functions Ani{r) and Bni{r) have rather 
complicated expressions: 



A^Ar) 



B„Ar) 



a 



nl 



nl 



^^Jii<l'nir)-l{l + l)KjJ^+ 
Qni dr Kir 



+ /^n.^^l/z(g„V)-/(/ + l)/^.z^ 



kf^ir dr 



+ n 



nl 



y iiQnir) 
Qnir 



(4.59) 



(4.60) 



Here k^iR and gf^i? are dimensionless eigenvalues , and they are the solution to a 
rather complicated algebraic equation for the frequencies u = — see for details. 
In and (|;60D I have set 



nl 



Ctq'ni 



D 



'^nl 



E 



D„i = 



(4.61) 



^ll^nl '^nl ^ll^nl 

and introduced the normalisation constant C„/, which is fixed by the orthogonality prop- 
erties 



/ i'^n'l'm'Y ■ i'^nlrn) Qod^X = M ^nn'^W^n 



(4.62) 



where M is the mass of the hollow sphere: 
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Equation ( ^.62| ) fixes the value of Cni through the radial integral 



R 



l{l + l)B. 



An 

r'^dr = —go (1 - ^^)R^ 



(4.63) 



(4.64) 



as can be easily verified by using well known properties of angular momentum operators 
and spherical harmonics. 1 shall later specify the values of the different parameters 
appearing in the above expressions as required in each particular case which will in 
due course be considered. As seen in reference a scalar-tensor theory of GWs such 
as JBD predicts the excitation of the sphere's monopole modes as well as the m = 
quadrupole modes. In order to calculate the energy absorbed by the detector according 
to that theory it is necessary to calculate the energy deposited by the wave in those 
modes, and this in turn requires that 1 solve the elasticity equation with the GW driving 
term included in its right hand side. The result of such calculation was presented in full 
generality in reference 0, and is directly applicable here because the structure of the 
oscillation eigenmodes of a hollow sphere is equal to that of the massive sphere — only 
the explicit form of the wavefunctions needs to be changed. I thus have 



m=—l 



(4.65) 



where G'^^'^\uJni) is the Fourier amplitude of the corresponding incoming GW mode, and 



JnO 



Jn2 



' M 



R 



(4.66) 
(4.67) 



for monopole and quadrupole modes, respectively, and Ani{r) and Bniir) are given by 
( |4.59| ). Explicit calculation yields 



'r 

bn2 

R 



3 a 



nO 



47r 1 - ^3 
3 Cn2 

47r 1 - ^3 



A{R) - c:3A(a) 
S(i?) - ^3s(a) 



(4.68) 
(4.69) 



with 



A{z) ^ (4.70) 

QnoR QnoR 

m - ''^'"-f - 3i^.. 'f^'f + Bra - 3£)„. ^ (4.71) 

fe-R kn2R Qn2R kn2R 

The absorption cross section, defined as the ratio of the absorbed energy to the 
incoming flux, can be calculated thanks to an optical theorem, as proved e.g. by Weinberg 
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28|. According to that theorem, the absorption cross section for a signal of frequency 



uj close to (-J AT, say, the frequency of the detector mode excited by the incoming GW, is 
given by the expression 

a(u}) = 1- — — — 4.72 

where F is the linewitdh of the mode — which can be arbitrarily small, as assumed in the 
previous section — , and rj is the dimensionless ratio 

., = %^ = i^ (4.73) 

J- J- -C/osc 

where Pqw is the energy re- emitted by the detector in the form of GWs as a consequence 
of its being set to oscillate by the incoming signal. In the following I will only consider 

the case PqW = -Pscalar-tensor with [i, g 



scalar— tensor 



\Qkkiuj)\' + ^Q^uj)Q,,iuj] 



5c5 {2LUBD + 3) 

where Qij{uj) is the quadrupole moment of the hollow sphere: 



(4.74) 



Qijiuj) = I XiXj g{'K,uj) d x (4-75) 

J Antenna 

and ujbd is Brans-Dicke's parameter. 



5. Scalar-tensor cross sections 



Explicit calculation shows that -Pscaiar-tensor is made up of two contributions: 

-Pscalar— tensor Poo + P20 (5.76) 

where Pqo is the scalar, or monopole contribution to the emitted power, while P20 comes 
from the central quadrupole mode which, as discussed in [§] and 0, is excited together 
with monopole in JBD theory. One must however recall that monopole and quadrupole 
modes of the sphere happen at different frequencies, so that cross sections for them only 
make sense if defined separately. More precisely. 



CTnOK.^) 



107rr7„2 



r2 



/4 



[UJ 



(5.77) 
(5.78) 



where rjno and ?7„2 are defined like in ([4.73|) , with all terms referring to the corresponding 
modes. After some algebra one finds that 



{ubd + 2) (u; - ujnor + r„o/4 

"'^^^"^ - "^^(-.^. + 2)03 (.-..„2)^ + ry4 ('-'^^ 
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Here, I have defined the dimensionless quantities 



4vr2 

Hn = TTJ— -AK^hnof (5.81) 

9 (1 + CTp) 



87r 



2 



15 l + (Tp 



(A;„2&n2)' (5.82) 



where cxp represents the sphere material's Poisson ratio (most often very close to a value 
of 1/3), and the h^i are defined in ([4.68|) ; vs is the speed of sound in the material of the 
sphere. 

In tables |l| and ^ I give a few numerical values of the above cross section coefficients. 



Table 1. Eigenvalues k^^R, relative weights Dno and iJ„ coefficients for a hollow sphere 
with Poisson ratio ap=l/3. Values are given for a few different thickness parameters 





n 




DnO 




0.01 


1 


5.48738 


-.000143328 


0.90929 




1 


12.2332 


-.001.59636 


0.14194 




2 


18.6321 


-.00558961 


0.05926 




4 


24.9693 


-.001279 


0.03267 


0.10 


1 


5.45410 


-0.014218 


0.89530 




1 


11.9241 


-0.151377 


0.15048 




2 


17.7277 


-0.479543 


0.04922 




4 


23.5416 


-0.859885 


0.04311 


0.15 


1 


5.37709 


-0.045574 


0.86076 




2 


11.3879 


-0.434591 


0.17646 




3 


17.105 


-0.939629 


0.05674 




4 


23.605 


-0.806574 


0.05396 


0.25 


1 


5.04842 


-0.179999 


0.73727 




2 


10.6515 


-0.960417 


0.30532 




3 


17.8193 


-0.425087 


0.04275 




4 


25.8063 


0.440100 


0.06347 


0.50 


1 


3.96914 


-0.631169 


0.49429 




2 


13.2369 


0.531684 


0.58140 




3 


25.4531 


0.245321 


0.01728 




4 


37.9129 


0.161117 


0.07192 


0.75 


1 


3.26524 


-0.901244 


0.43070 




2 


25.3468 


0.188845 


0.66284 




3 


50.3718 


0.093173 


0.00341 




4 


75.469 


0.061981 


0.07480 


0.90 


1 


2.98141 


-0.963552 


0.42043 




2 


62.9027 


0.067342 


0.67689 




3 


125.699 


0.033573 


0.00047 




4 


188.519 


0.022334 


0.07538 



As already stressed in reference [27|, one of the main advantages of a hollow sphere 



is that it enables to reach good sensitivities at lower frequencies than a solid sphere. 



Table 2. Eigenvalues k^2^y relative weights Kn2, Dn2, Dn2 and F„ coefficients for a 
hollow sphere with Poisson ratio ap = 1/3. Values are given for a few different thickness 
parameters q. 





n 




Kn2 


Dn2 


Dn2 




Fn 


0.10 


1 


2.63836 


0.855799 


0.000395 


-0.003142 


2 


.94602 




2 


5.07358 


0.751837 


0.002351 


-0.018451 


1, 


.16934 




3 


10.96090 


0.476073 


0.009821 


-0.071685 


0, 


.02207 


0.15 


1 


2.61161 


0.796019 


0.001174 


-0.009288 


2, 


.86913 




2 


5.02815 


0.723984 


0.007028 


-0.053849 


1, 


.24153 




3 


8.25809 


-2.010150 


-0.094986 


0.672786 


0, 


.08113 


0.25 


1 


2.49122 


0.606536 


0.003210 


-0.02494 


2, 


.55218 




2 


4.91223 


0.647204 


0.019483 


-0.13867 


1, 


,55022 




3 


8.24282 


-1.984426 


-0.126671 


0.67506 


0, 


.05325 




4 


10.97725 


0.432548 


-0.012194 


0.02236 


0, 


.03503 


0.50 


1 


1.94340 


0.300212 


0.003041 


-0.02268 


1, 


.61978 




2 


5 06453 


745258 


0051 33 


-0 02889 


2 


.29572 




3 


10.11189 


1.795862 


-1.697480 


2.98276 


0, 


.19707 




4 


15.91970 


-1.632550 


-1.965780 


-0.30953 


0, 


.17108 


0.75 


1 


1.44965 


0.225040 


0.001376 


-0.01017 


1, 


.15291 




2 


5.21599 


0.910998 


-0.197532 


0.40944 


1, 


,82276 




3 


13.93290 


0.243382 


0.748219 


-3.20130 


1, 


,08952 




4 


23.76319 


0.550278 


-0.230203 


-0.81767 


0, 


,08114 


0.90 


1 


1.26565 


0.213082 


0.001019 


-0.00755 


1, 


,03864 




2 


4.97703 


0.939420 


-0.323067 


0.52279 


1, 


,54106 




3 


31.86429 


6.012680 


-0.259533 


4.05274 


1, 


,46486 




4 


61.29948 


0.205362 


-0.673148 


-1.04369 


0, 


,13470 



For example, a hollow sphere of the same material and mass as a solid one (<j = 0) has 
eigenfrequencies which are smaller by 

o^niW =w„i(^ = 0) (1-^3)^/3 (5.83) 

for any mode indices n and /. I now consider the detectability of JBD GW waves coming 
from several interesting sources with a hollow sphere. 

The values of the coefficients Fn and Hn, together with the expressions ( p.77| ) for 
the cross sections of the hollow sphere, can be used to estimate the maximum distances 
at which a coalescing compact binary system and a gravitational collapse event can be 
seen with such detector. I consider these in turn. 

By taking as a source of GWs a binary system formed by two neutron stars, each 
of them with a mass of mi =m2 = 1.4 M©. The chirp mass corresponding to this sys- 
tem is Mc = (mim2)^^^ {nii + 7712)"^^^ = 1.22 Mq, and z/[5 cycles] = 1270 Hz. Repeating the 
analysis carried on in Section three I find a formula for the minimum distance at which a 
measurement can be performed given a certain signal to noise ratio (SNR), for a quantum 
limited detector 
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r[UJnOj 



5 -21/3 



32 (fiBD + 2)(12^]BB + 19) 

5 -21/3 1 G5/3jVf5/3 

192 (^]Bi5 + 2)(12^]BZ) + 19) 



nl/2 



huti^SNR 



huSSNR 



1/2 



(5.84) 
(5.85) 



For a CuAL sphere, the speed of sound is f5 = 4700 m/sec. I report in table 
^ the maximum distances at which a JBD binary can be seen with a 100 ton hollow 
spherical detector, including the size of the sphere (diameter and thickness factor) for 
SNR = 1. The Brans-Dicke parameter Qbd has been given a value of 600. This high 
value has as a consequence that only relatively nearby binaries can be scrutinised by 
means of their scalar radiation of GWs. A slight improvement in sensitivity is appreciated 
as the diameter increases in a fixed mass detector. Vacancies in the tables mean the 
corresponding frequencies are higher than z/p cycles]- 



Table 3. Eigenfrequencies, sizes and distances at which coalescing binaries can be seen 
by monitoring of their emitted JBD GWs. Figures correspond to a 100 ton CuAl hollow 
sphere. 



<; 


$ (m) 


1^10 (Hz) 


1/12 (Hz) 


r(i^io) (kpc) 


r(j/i2) (kpc) 


0.00 


2.94 


1655 


807 




29.8 


0.25 


2.96 


1562 


771 




30.3 


0.50 


3.08 


1180 


578 


55 


31.1 


0.75 


3.5 


845 


375 


64 


33 


0.90 


4.5 


600 


254 


80 


40 



The signal associated to a gravitational collapse can be modeled, within JBD theory. 



as a short pulse of amphtude b, whose value can be estimated as [29 



10 



-23 



500 



n 



MrJ \ 



f 10 Mpc\^ 



(5.86) 



Table 4. Eigenfrequencies, sizes and distances at which coalescing binaries can be seen 
by monitoring of their emitted JBD GWs. Figures correspond to a 3 metres external 
diameter CuAl hollow sphere. 





M (ton) 


i^io(Hz) 


1/12 (Hz) 


r{iyio) (kpc) 


r{iyi2) (kpc) 


0.00 


105 


1653 


804 




33 


0.25 


103.4 


1541 


760 




31 


0.50 


92 


1212 


593 


52 


27.6 


0.75 


60.7 


997 


442 


44.8 


23 


0.90 


28.4 


910 


386 


32 


16.3 
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where is the collapsing mass. 

The minimum value of the Fourier transform of the amplitude of the scalar wave, 
for a quantum limited detector at unit signal-to-noise ratio, is given by 

where K^, = 2if„ for the mode with / = and Kn = Fn/3 for the mode with Z = 2, m = 0. 

The duration of the impulse, r ^ l//c, is much shorter than the decay time of the 
nl mode, so that the relationship between b and b{ujni) is 

b ^ \b{ujni) \ fc atfrequency Uni = 27r/c (5.88) 
so that the minimum scalar wave amplitude detectable is 

Let us now consider a hollow sphere made of molibdenum, for which the speed of 
sound is as high as fs = 5600 m/sec. For a given detector mass and diameter, equation 
(|5.89| ) tells us which is the minimum signal detectable with such detector. For exam- 
ple, a solid sphere of M = 31 tons and 1.8 metres in diameter, is sensitive down to 
^min = 1-5 ■ 10~^^. Equation ( 5.86| ) can then be inverted to find which is the maximum 



distance at which the source can be identified by the scalar waves it emits. Taking a 
reasonable value of Qbd = 600, one finds that r(z/io) ~0.6 Mpc. 

Like before, I report in tables ||, ^ and |^ the sensitivities of the detector and 
consequent maximum distance at which the source appears visible to the device for 
various values of the thickness parameter ^. In table ^ a detector of mass of 31 tons has 
been assumed for all thicknesses, and in tables ||, |^ a constant outer diameter of 3 and 
1.8 metres has been assumed in all cases. 



Table 5. Eigenfrequencies, maximum sensitivities and distances at which a gravitational 
collapse can be seen by monitoring the scalar GWs it emits. Figures correspond to a 31 
ton Mb hollow sphere. 



<; 


(j> (m) 


1/10 (Hz) 


|6|min (10-22) 


r(zyio) (Mpc) 


0.00 


1.80 


3338 


1.5 


0.6 


0.25 


1.82 


3027 


1.65 


0.5 


0.50 


1.88 


2304 


1.79 


0.46 


0.75 


2.16 


1650 


1.63 


0.51 


0.90 


2.78 


1170 


1.39 


0.6 
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Table 6. Eigenfrequencies, maximum sensitivities and distances at which a gravitational 
collapse can be seen by monitoring the scalar GWs it emits. Figures correspond to a 1.8 
metres outer diameter Mb hollow sphere. 





M (ton) 


i^io (Hz) 


|6|min (10-22) 


r{vw) (Mpc) 


0.00 


31.0 


3338 


1.5 


0.6 


0.25 


30.52 


3062 


1.71 


0.48 


0.50 


27.12 


2407 


1.95 


0.42 


0.75 


17.92 


1980 


2.34 


0.36 


0.90 


8.4 


1808 


3.31 


0.24 
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